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Abstract 



In the paper, we consider a special coupled forward-backward stochastic differential system (FB- 
SDS) which is associated to the viscous incompressible Navier-Stokes equation and provides a prob- 
abilistic solution to the latter via the Feynman-Kac formula. With a probabilistic method, we first 
£N| . prove the existence and uniqueness of the solution to the FBSDS. Then under the same conditions, 

we verify that this solution leads to a unique local strong solution to the associated Navier-Stokes 
equation, and for both cases of the small Reynolds number and dimension two, we further give the 
O |' global strong solutions. Our probabilistic representation formula for the solution to the Navier-Stokes 

equation involves neither integral nor gradient operators, and it can serve to generate a Monte Carlo 
solution to the viscous incompressible Navier-Stokes equation. 



Keywords: forward-backward stochastic differential system, Navier-Stokes equation, Feynman-Kac 
formula. 



: 1 Introduction 

CN , 

The standard deterministic Navier-Stokes equation describes the evolution of the velocity field of an 
incompressible, viscous fluid moving in a domain of R d (d = 2 or 3 throughout this work) , and takes the 
C*~) . following form: 

Si ( d t u-^Au+ (u-V)u + Vp + f = 0, r>0; 

V • u = 0, w(0) = mo, 

where u is the d-dimensional velocity field of a fluid, p is the pressure field, v £ (0,oo) is the viscosity 
coefficient, and / is the external force which, without any loss of generality, is taken to be divergence 
free. Let T £ (0,oo) be a real sufficiently big number. If (u,p) solves the initial Cauchy problem 
then (u, p) defined by the following time-reversing transformation 

u(t,x) = —u(T — p(t,x) = p(T — t,x), f(t,x) = f(T — t,x), for t < T, 

solves the following terminal Cauchy problem: 

d t u+%Au + (u-V)u + Vp + f = 0, t<T; 

V-u = 0, u(T) = G:= -it 0) [ ' ' 

which is also called (backward) Navier-Stokes equation due to its equivalence to the former. 
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It is extremely hard to completely solve the Navier- Stokes equation (jl.2p or and it still remains to 

be a striking important open problem. Charles Fefferman in his well-known article [17] finally commented, 
"Standard methods from PDE appear inadequate to settle the problem. Instead, we probably need some 
deep, new ideas." In the paper, we develop a probabilistic methodology to study the Navier-Stokcs 
equation. More precisely, the Navier-Stokes equation (11.21) is associated to the following coupled forward- 
backward stochastic differential system (FBSDS): 

dX s (t,x) =Y s (t,x)ds + V^dW s , se [t,T]; 
X t (t,x) = x; 

-dY s {t,x) = [f(8,X,(t,x))+Y (8,X a (t,x)j\ ds-V^Z s (t,x)dW s , s€ [t,T]; 

Y T (t,x)=G(X T (t,x)); 

d ^ 1 

-dY s (t, x) = £ ^Yt l Yt(t, * + B s ) (BL - B* ) [B{ - BL) B § ds 

»>i=i 

- dM s , s e (0,oo); 
Y oo (t,x)=0. 

Here and in the following, B and W are two independent cZ-dimensional standard Brownian motions. 
The drift part of {Y s (t,x),s € [i,? 1 ]} (see the third equality of FBSDS (11.31) ) at time s depends on 
Y , and that of {Y s (t,x),s £ [t,T]} depends on Y t (t,x + B s ), which make our system (ll.3[) strikingly 
different from the conventional coupled forward-backward stochastic differential equations (FBSDEs) 
(see [TJ HD1 H3 E2 EZl EH S3). Furthermore, both backward stochastic differential equations (BSDEs) 
in FBSDS (|1.3j) are defined on two different time-horizons [t,T] and (0,oo). The i? m -solution (to > d/2, 
see Definition I3.2[) (X,Y, Z,Y Q ) to FBSDS (jl.3|) will be connected to the strong solution (u,p) to the 
Navier-Stokes equation (|1.2| in the following manner: 



Y 3 (t,x) = u(s,X s (t,x)), Z s {t,x) = Vu(s,X s (t,x)), and Y (t,x) = Vp(t,x), (s,x) £ [t,T] x R d . 

FBSDEs have been connected to a system of semi-linear parabolic partial differential equations (PDEs) 
(see among many others [TJ[3ni[2iB[2Ill251II3])- However, Navier-Stokes equation (|1.2[) usually goes beyond 
that context, as it has the nonlocal constraint V • u = 0. Our difficulty is two-fold: one is nonlinearity 
and the other is the nonlocal constraint. To attack the divergence- free constraint, we introduce the last 
BSDE in the infinite time interval. To concentrate our attention on the primary connection between 
FBSDS (|1.3|) and Navier-Stokes equation ()1.2|) . FBSDS (|1.3p is discussed from an analytic point of view 
only in a Hilbert space, though it might be addressed in other spaces like Holder space. 

There is a long history on the formalisms to represent solutions of PDEs as the expected functionals of 
stochastic processes. We only mention here those concerning a deterministic incompressible Navier-Stokes 
equation. The velocity field is related to the vorticity field in a linear fashion by the Biot-Savart law, 
and moreover, the analysis of the vorticity field is fundamental to the issues like the possible emergence 
of singularities (for instance, see [3J |3D]). For the two-dimensional case, since the vorticity obeys a 
Fokker-Planck type parabolic PDE, the random vortex method was formulated by Chorin [5] who used 
random walks and a particle limit to represent the vorticity field, and Busnello [6] used the Girsanov 
transformation to give a probabilistic representation of the vorticity field. The latter work was further 
extended by Busnello, Flandoli and Romito [7] to the three-dimensional case, where the vorticity field 
turns out to satisfy a parabolic PDE with an additional stretching term. As noted by Busnello, Flandoli 
and Romito , a partially similar representation formula for the vorticity field of the three-dimensional 
Navier-Stokes equations had been given before by Esposito et al. [TBI E] but without the probabilistic 
representation for the Biot-Savart law. Note that a probabilistic interpretation for the Biot-Savart law 
was given by Busnello, Flandoli and Romito [7] and Busnello [6J, where the Bismut-Elworthy-Li formula 
is used so that the velocity can be recovered from the vorticity through probabilistic approaches. Le 
Jan and Sznitman [24 interpreted the Fourier transformation of the Laplacian of the three-dimensional 



2 



velocity field in terms of a backward branching process and a composition rule along the associated tree, 
and got a new existence theorem, and their approach was extensively studied and generalized by others 
(see, for instance [4] 132]). 

Recently, Constantin and Iyer [£l HQ] and Iyer [22 [221 E3] derived a stochastic representation for the 
incompressible Navier-Stokes equations based on stochastic Lagrangian paths and gave a self-contained 
proof of the existence. Later, Zhang [44] considered a backward analogue and provided short elegant 
proofs for the classical existence results. However, their probabilistic representations were not so com- 
plete as those of [3 [HI [Ml 122] in the sense that they retain the Leray-Hodge projection P and the 
gradient operator (see [HI [33 [23 1221 1221 EI] ) • Ours involves neither integral nor gradient operators, and 
it can serve to generate a Monte Carlo solution to the Navier-Stokes equation. We also mention that 
Cruzeiro and Shamarova [11] established a connection between the strong solution to the spatially peri- 
odic Navier-Stokes equations and a solution to a system of FBSDEs on the group of volume-preserving 
diffeomorphisms of a flat torus, and that Qiu, Tang and You [38] considered a similar non-Markovian 
FBSDS to ours (|1.3p in the two-dimensional spatially periodic case, and studied the well-posedness of 
the corresponding backward stochastic PDEs. The list of literature on probabilistic approaches to the 
Navier-Stokes equations is not exhausted here, and there are many others. 

The rest of this paper is organized as follows. In Section 2, we introduce notations and functional 
spaces, and recall auxiliary results. In Section 3, a lemma on equivalent norms and the definition of the 
solution to FBSDS (ll.3[) are given first and then the FBSDS is connected to the Navier-Stokes Equation, 
which constitutes our main result of Theorem 13.31 In Section 4, we prove Theorem 13.31 In Section 5, 
we give global existence of the solution for the small Reynolds number and the two-dimensional cases, 
respectively. Finally in Section 6 as an appendix, we prove Lemmas 13.11 and 14.11 



2 Preliminaries 

Let (Q, {^t}t>0j P) be a complete filtered probability space on which are defined two d-dimensional 
standard Brownian motion W — {Wt : t £ [0, oo)} and B — {B t : t £ [0, oo)} such that {J^t}t>o is 
the natural filtration generated by W and B, and augmented by all the P-null sets in # . By {J^}t>o 
and {^ B }t>o, we denote the natural filtration generated by W and B respectively, and they are both 
augmented by all the P-null sets. 3? is the a- Algebra of the predictable sets on SI x [0, T] associated with 
{&t}t>o- 

Denote by | • | (respectively, (•, •) or •) the norm (respectively, scalar product) in finite-dimensional 
Hilbert space such as R, M fe , R fex ' where k,l are positive integers and 



1/2 



1/2 



and \y\ := [ £ £; 

;=i j=i 



for (x, y) £ 



i k x R kxl . 



For each Banach space (X, || • \\x) and real q £ [1, oo], we denote by S q ([t, r]; X) the set of ^-valued, 
^-adapted and cadlag processes {X s } s£ [ t ^ T ] such that 



l^nU«([t,T] ; ;v) : - E 



sup || X 



8\\X 



1/'/ 



< OO. 



L q ^(t,T; X) denotes the set of (equivalent classes of) AWalued predictable processes {A^ s } se [ tiT ] such that 

1/9 



\x\ 



E 



\X s f x ds 



< oo. 



Both (S 9 ([t, t];X), II ■ \\ S i([t,T]x)) and f L%(t, r; X), \\ ■ \\ L ^ t , T -x)) are Banach spaces. 
Define the set of multi-indices 



A := {a — (ai, • • • , ad) ■ ol\, ■ ■ ■ ,ad are nonnegative integers}. 
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For any a e A and x = (xi, • • • , xj) e R d , denote 

d 



dx^dx^ 2 ■■■dx° d ' 

For differentiable transformations 0,-0 on R d , define the Jacobi matrix V</> of </>: 

d Xl 4?,d X2 (f> 2 , ■ ■ ■ ,d Xd 4> 2 



V0 



(2-1) 



V d Xl 4> d ,d X2 <p d ,--- ,d Xd <f> d J 



whose transpose is denoted by V*</>, the divergence div</> = V • 0, and the matrix 



(2.2) 



Now we extend several spaces of real- valued functions to those of vector- valued functions. For a 
positive integer number I and fc, we denote by (R* ; R k ) (respectively, C^ o (0; M fe ) for each open set 
O C R ) the set of all infinitely differentiable IR^-valued functions with compact supports on M 1 (O, 
respectively) and by £)'(R l ;R k ) the totality of all the Revalued general functions with each component 
being Schwartz distribution. For simplicity, we write and $' for the case I = k = d. On R d we denote 
by (S*", respectively) the set of all the Revalued functions whose elements are Schwartz functions 
(tempered distributions, respectively). We shall denote by (•, •) not only the duality between and S' 
but also the duality between 5? and Sf' . Then the Fourier transform T{f) of / € is given by 

T(f)(0 = (2n)- d / 2 f cxp (-V=l(x, 0)f(x) dx, £ € R d , 

JR d 

and the inverse Fourier transform F~ 1 {f) is given by 

F-\f){x) = {2n)- d ' 2 f exp (V=T(x, £>)/(0 ^, x E R d . 

Extended to the general function space y, the Fourier transform defines an isomorphism from onto 
itself. As usual, for each s e R and f e =5"', we denote the Bessel potential I s (f) := (1 — A) s / 2 / = 
J--i((l + |£H*/2.F (/)(0) . 

For each positive integer Z, 1 < g < oo and m = 0, 1, . . . by L q (R l ) and H m ' q (R l ) (L q and H m < q , with 
a little notional abuse), we denote the usual R'-valued Lebesgue and Sobolev spaces on R d , respectively. 
H m ' q is equipped with the norm: 



= { (\\4>\\l« + Y^\=AD a 4>\\ q Lq ) 

ll0lk~ + EM=iP < VlU~, 



1/9 



which is equivalent to the norm: 



||(l-A)T0|| i4) 0eif ro ' 9 , forge (l,oo). 



Both norms will not be distinguished unless there is a confusion. Furthermore, using the Bessel potentials, 
we define the Soblev space H m ' q := I- m (L q ) for m e Z \ (0 U N) and g € (1, oo). In particular, for the 
case of g = 2, iJ m - 2 is a Hilbcrt space with the inner product: 

(0, V)m := / (I m /2(t>{x), I m/2 ip{x))dx, <j>, V € i? m ' 2 . 
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We define the duality between H s ' q and H r,q for q G (1, oo) and q' = q/(q— 1) as: 
(0, 1>).,r := / (/./a0(*), I r/ a^(x)) di, G i/ s ' 9 > G W. 



For simplicity, we write the space H m and the norm || • || m for H m ' 2 and || • \\ m ,2, respectively. 
Define 

9 a := {0 G C c °° : V ■ - 0} . 

Denote by iJ™' 9 the completion of ^ under the norm || • || m , g , which is a Banach subspace of H m ' q . 

Now we introduce several spaces of continuous functions. For each positive integer /, nonnegative 
integer k and domain O C R d , we denote by C(0,R l ), C k {0,R l ) and C k > s (0, R l ) with S G (0,1) the 
continuous function spaces equipped with the following norms respectively: 

fe 

H0llc(O,E') : = SU P 100*01, H0llc*=(O,K') : = H0llc(O,R ! ) + Z] ll^^llcCO.R')' 

n ,,, || j II . \- o 

\\ ( P\\c k -> s (0,R 1 ) ■= H0llc*(0,R ! ) + 2^ SU P 1 _ [g J 

with the convention that C (O,R') = C(0,R l ). Whenever there is no confusion, we write C(R d ),C k , 
and C M for C(R d ,R l ),C k (R d , R' ) and C k ' 5 {R d ,R l ), respectively. We define C°°(R d ) := n™ =1 C k (R d ). 

In an obvious way, we define spaces of Banach space valued functions such as C(0, T; H m ' q ) and 
L r (0, T; H m ' q ) for m £ Z,r,q G [1, oo], and related local spaces like the following ones: 

L r loc (T ,T;H m - q ):= (J L r (Ti,T; H m ' q ), Ci oc ((T ,T]; H m ' q ) :— (J C([Xi, T]; iF™' 9 ). 

Tie(T ,T] Tie(T ,T] 

We have the following properties on Sobolev spaces, whose proof is omitted. 

Lemma 2.1. There holds the following assertions: 

(i) the space H n , n > d/2 + k, k G Z + U 0, is continuously embedded into the space C k ' S for any 
5 G (0, (n — d/2 — fc) A 1), i.e., i/iere exists a constant C > sucft i/iai 

||0|| c m <C||0||„, V0 e if"; 

(ii) if 1 < r < s < oo and m,n G Z satisfying - — m = - — n, then H n ' r is continuously embedded 
into H m,s , i.e., there exists a constant C > suc/i t/iat 

||0||m,. < C||0|U,r, V0GiT 1 ' r ; 

(mj /or m G Z + U {0}, £/iere exists a constant C > swc/i £/iai, /or any 0, ip G n i? m , 

IWIU<C{||0||L-.||i? m ^||L» + 11^11^11^011^}, 

£ ||C Q (0^)- 00^11^ <C{l|V0|| L ^|| J D m -VllL 2 + ll^ m 0llLHIV'IU~}, m>L; 

0< |a| < rn 

(iv) for any s > d/2, H s is a Banach algebra, i.e., there exists a constant C > such that, 

\\<H>\\. < CUW.Uh, V0,^ei/ s . 



The first two assertions of Lemma 12.11 are borrowed from the well-known embedding theorem in 
Sobolev space theory (see [H]) and the others are referred to [301 Lemma 3.4, Page 98]. For simplicity, 
we shall denote by <—> the embedding relationship, i.e., by A S- B we mean that normed space (A, || • \\a) 
is embedded into (£?, || • \\b) with a constant C such that 

||/||b<C||/|U, v/g A 
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Remark 2.1. Note that d — 2 or 3 throughout this work. For any h,g £ H 2 , we have 

\\h-g\\ 2 

= \\h- g\\l + J2\\ d *> h - 9\\o + it\\ h - d ^o 

i=l i=l 

<^{||/ l ||| ?(H , ) ||ff||§ + ||Vft||g''||Vfe||?- !W, ||fl||^||5||?- V + MINIMI?} 
( using Gagliardo-Nirenberg Inequality (see [T8l l28l I3T] )) 

<c||ft||lNI? s 

where /3 := 1 — d/4, and -ff 2 is embedded into C°' S for some <5 £ (0, 1). In view of Lemma [2 .11 we have 
for any integer m > d/2, 

INIL-l < C("*,d)IWUIslL-i> Vft £H m ,g£ H m -\ 
Lemma 2.2. There hold the following : 

K(0- V)V, V-)m-i, m+ i| < C|MUlVllm. m>2,<j>£H™^£ h: i+1 . (2.3) 

and 

K^-v^.^l.rciiv^lU-^W^ m>2,^c +1 - (2-4) 

Proof. The first inequality (|2.3p is referred to [26 . For the reader's convenience, we give a simple proof 
of the second inequality (|2.4[) . 

For any multi- index a such that \a\ < m, we have 

£ Q ((0-V)</>) = (0-V)i} Q + ]T C(a,/3)(D^-V)D a -^. 

0</3<q 

Since ((</> ■ V)D a <fi, _D Q 0)o = 0, we have from Assertion (iii) of Lemma \2. II that 

0<|a|<m 0<^<a 
d 

2 



< C(m, d)||V0|| L oc||0|| 



□ 



3 Connection between FBSDS and the Navier- Stokes Equation 

3.1 FBSDEs with coefficients in Sobolev Spaces 

Assume that v > and that 

b £ C([0, T]-H m ) n L 2 (0, T; H m+1 ), 4> £ L\0, T; L 2 {R d )), i/> £ L 2 {R d ), (3.1) 

for some integer m > d/2. Consider the following FBSDE: 

' dX s (t,x) = b(s,X s (t,x))ds + y/vdW s , s £ [t,T]; 
XAt, x) = x\ 

r- (3-2) 

-dY s (t,x)=<t>(s,X s (t,x))ds-VvZ s (t,x)dW s , s£ [t,T]; 
Y T (t,x)=iP(X T (t,x)). 



Since H m <-> C ' 5 and H m+1 <-> C 1 ' 5 for m > d/2, in view of [27J Theorems 3.4.1 and 4.5.1], the forward 
SDE is well posed for each (t,x) £ [Tq,T] x R d , and the unique solution in relevance to the initial data 
(t, x) £ [To, T] x R d defines a stochastic flow of homeomorphisms. Since the function (j> is only measurable, 
the following lemma serves to give a clear meaning of the composition <^>(s, X s (t, x)). 

Lemma 3.1. Assume that m > d/2 andb £ C([0, T]; H m ) n L 2 (0, T; H m+1 ). Then there are two positive 
constants n and K which only depend on \\ div b\\ z,i(o,T;Z,°°)> such that for all t £ [0, T], s £ [t, T], (</?, 77) £ 
L 1 (R') x L 2 ([0,T] x R d ; R'), and Z e Z+ ; we have 

K|Mlii(R')</ S[k(^.(t,as))l]<to<«'lkll£i(R«) ) (3-3) 

^IhlliH^TlxR') < / / ^[Ms^&a))!] dsda; < -K"||7?||ii([ t) T]xR')- (3- 4 ) 



Our preceding lemma weakens the assumptions on b of [2J Theorem 14.3], where b(t, ■) = b(-) is 
time invariant and is required to lie in the more regular space C 1 (R d ) n W 1,oa (R d ). Since b(t,x) is not 
necessarily uniformly Lipschitz continuous in x, the stability of X with respect to the coefficient b has 
to be proved very carefully and the proof of [2] has to be generalized accordingly, we give a probabilistic 
proof in the appendix. 

Remark 3.1. From Lemma [6. 11 we see that Lebesgue's measure transported by the flow {X s (t, x),s £ 
[t, T]} results in a group of measures {/i s , s £ [t, T]} satisfying for any Borel measurable set A C R d , 



H.(A) = / E[l A (X s (t,x))} dx. 

These measures are all equivalent to Lebesgue measure and the exponential rate of compression or dilation 
are governed by the divergence of b. This is similar to that of a system of ordinary differential equations 
(see [IS])- On the other hand, thanks to Lemma [3.11 our FBSDE (|3. 21) makes sense under assumption 
(|3~T|) . i.e., the forward SDE is well posed for each (t, x) £ [0, T] x R d and for each t £ [0, T] the BSDE is 
well posed for almost every x £ R d . 

3.2 Definition of the solution to system (11. 3p 

First, let us consider the following trivial BSDE on [0, 00]: 

-dS t = g t dt-dM t ; ^=0. (3.5) 

In the conventional sense (see [51 [T21 [TH [3S1 |37]), a solution of BSDE is always defined as a pair of 
processes (S, M) and the process M serves to guarantee the adaptedness of S. However, since FBSDS 
(|1.3[) only involves Yq rather than the whole process Y, we are interested only in determining So instead 
of the whole process S for such kind of BSDEs on [0, 00]. 

Different from [5J [T21 [TH [3H1 137] , we assume that {gt, t £ (0, 00)} is an Revalued J^-adapted process 
such that 

/■OO poo 

/ E \\g t \] dt < 00, Ve > and lim / E [g t ] dt < 00. (3.6) 

The class £ is defined as the totality of BSDEs (|3 . 5|) with the drift satisfying the preceding assumption. 
Each BSDE (|3.5|) lying in £ can be written into the integral form: 

/•OO 

S t = I g s ds- M t , 



where as we can not solve it in the whole time interval [0,oo) at once, we solve it with a pair (S,M) 
instead of (S,M), with E[M t \^ t B ] = 0, Vi £ (0,oo). Therefore, we have 



S t = E 



/oo poo 
g s ds\^ t B \ , M t = J g s ds - S u Vt £ (0, 00). 
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Sq is defined as follows 



Sq := \imE [S e ] , 



(3.7) 



while Mq is not necessary for us to know its meaning in this work. 

Now, our trivial BSDE p.5[) makes sense and the definition of the solution does not conflict with 
the common cases (see [5] [12, 14, 36, 37 ). Further, we may consider more general nonlinear cases such 
as g s — g(s, S s , M s ) which may be used to describe more general operators like in Lemma 13.21 below. 
However, we do not seek such a generality in this paper. 

Definition 3.1. We say that (X,Y, Z,Yq) is a solution (local solution, respectively) of FBSDS (|1.3|) if 

for almost every (t,x) £ [0,T] x R d ((t,x) £ (T ,T] x R d for some T Q £ (0,T), respectively), the BSDE 
on the infinite time interval belongs to class £, and 

(X.(t,x),Y.(t,x),Z.(t,x)) £ S 2 ([t, T];R d ) x S 2 ([t, T]:R d ) x L%(t,T;R dxd ) 

such that the first two stochastic differential equations of (jl.3l) on [0, T] (any subinterval [Xi,T] with 
Ti £ (Tq,T), respectively) hold almost surely. 

If both the BSDE on the infinite time interval [0, oo] and its unknown variable Y do not appear in 
FBSDS (| 1 .3(1 . then Definition 13.11 becomes automatically the definition of the solutions (local solutions, 
respectively) for an FBSDE by abandoning all the requirements on both Yq and the BSDE on the infinite 
time interval. 

To specify the regularity of the solutions, we introduce the following definition. 

Definition 3.2. (X, Y 1 Z, Yq) is called a (local, respectively) i? m -solution of FBSDS (|1.3p if it is a solution 
(local solution, respectively) of FBSDS (|1.3[) (on some time interval (T ,T], respectively) with Y £ 
L 2 (0,T;ii' m - 1 ) (Y Q £ Lf QC (T , T; if™" 1 ), respectively) and for each t £ [0,T] (t £ [T ,T], respectively) 
and almost every x £ R d , {Y s (t,x), s £ [t,T]} £ L°£(t, T; R d ). 

Remark 3.2. In FBSDS (|1.3|) . the regularity of Yo depends on that of Y t (t, x). However, Proposition ^. 31 
below shows that the regularity of Yo dominates that of (Y, Z). This explains why our ff m -solutions of 
FBSDS (|1.3[) only require the regularity of Y in Definition 13. 21 On the other hand, Definition 13 . 21 req uires 
the uniform boundedness of the unknown process Y(t, x). Indeed, in this work we only pursue the H m - 
solution (m > d/2, d = 2 or 3) of which Y.(t, x) is naturally uniformly bounded (see Theorem 13 . 31 below) ; 
moreover, the uniform boundedness is only required for the particular case m = 2 (see Proposition 14.31 
below). Compared with the standard assumptions for the solvability of FBSDEs (see [II [Tl [20l [29ll35l |43] 1 . 
it can be viewed as a compensation for the lack of continuity of the function f(s, x) with respect to x. 

3.3 Connections between FBSDS and Navier-Stokes equation 

Before we show the connections between the FBSDSs and Navier-Stokes equaitons, we give a probabilistic 
representation for an integral operator. 

Lemma 3.2. For <f>, ip £ H m with m > d/2 + I, define 

d 

i{x) := V(-A)- x dTO div{<j)®ij){x) = Vi-A^d^d^iftixWix)). 

Then, the following BSDE : 

-dY s (x) = ^^ j (x + B.) (El - BQ (bL - B|) S f ds 

- dM s , s £ (0,oo) ( 3 - 8 ) 

?oo(x) = 
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belongs to class £ for each x £ M. d and there holds 



C(x) =Y (x), Vxe 



(3.9) 



Proof. For m > d/2 + 1, H m is a Banach algebra embedded into if 2,7 for some 7 > d and also into 
C* u (R d ) for some S € (0, 1), so (/ftp € H m nH m ^ which implies d x id x j (ft (x)^ (x)) € if^M) nff ' 1 ^), 
i,j = l,--- , d. By the singular integral operator theory (see 39 ) and Lemma [2.11 we can check that 

£ e i? m_1 n c°- s (R d ). 

For each e > 0, we have 



E 



ds 



< c(U\\ m ,u\\ m ) 

COO -j^ 



27 

2s 3 



ds 



< C 



„3/2 



Thus, our BSDE (|3.8p is well-posed on the interval [e, 00] and we need only prove 



= limE YJx) 

eiO L . 



On the other hand, we have 



Y E {x) = E 



^^{x + B.) (BL - Bi) (Bi B± ) L% -/, 



,9? 



Careful calculations yield that, for each i,j,k= 1,2 
" 27 



E 

27 

: 2s 3 



_^(^ + S s )(5L-B|) (Bi -Bi)^ 



9 

2s2 



0V(a + y + z + r)yVr fe (27rs/3)~ 3<i/2 e 



(y)z j r k (2ns/3)- 3d/2 d y , 



3(i»r+ 



2 ^l~1 2 _i_l*.l 2 'i 



.|2 , 1 i2 , 1 i2 



dydzdr 



dydzdr 



9 

2s 2 



3 

2^ 
3 

2s" 



a yI (0 l ^)(2/)z J r fe (27rs/3)- 3d/2 e 
5 t/i (^^)(y)r fe (2 7 r S /3)- 3d / 2 e' 
^(^^')(y)r fc (2^/3)- 3d / 2 9 yJ < 



3(l«- 



-r-| 2 + |^| 2 + |r| 2 ) 



3|y- 



dydzdr 



■_ 3(|*| 2 + |i-| 2 ) 

" <9 zJ e 2s dydzdr 



2 = e 2s dzdydr 



3 

2s" 



^^ J (0V J )(y)r fc (47rs/3)' d/2 e~ 3 ' y '4T" 1 (27rs/3)~ d/2 e^^ riyrfr 



1 

2s" 



<9 yi 9^ V j ) (y) (4tts /3)' rf/2 e- 3 ' y '4r r| ( 27rs /3) - d / 2 d rfc e" 
fl„i 9 yJ (0 V ) (y) (47TS/3) ~ d/2 e- ^ (2tts/3) - d/2 d afc e~ ~ 
3 yi d yj (^^')(y)(27rs)- d / 2 ^ fe e-^ 
d x ,d x] (^)(y + x)y k {2^s)- d l 2 e-^ dy 



dydr 



drdy 



: —E [d x ,d x3 {^^){x + B s )B k ] , s > 0. 



(3.10) 



(3.11) 



(3.12) 



!) 



Moreover, we have 



\d xi d x3 {4> l ^){y + x)y k {2n S )- d / 2 t 



\dy 



< 



C 



1(1" 



(3.13) 



and 



<C\\d xi d xi {4>H j )\W q s-*-^, qe [1, 7 ] 

poo 

/ s^E [\d x ,d x3 {<t> l V){x + B s )B k s \] ds 

JQ 

\d xi d xi {<j> i ^){y + x)y k {2-Ks) 



-d/2 e -^ 



j dy ds 



= C ||9 x *S xi (^)l|o, 7 s-'-^ds + C / ||5^^(^)||o,i a 



(3.14) 



ds 



<c(\\d xi o xi (0 y ) 1 1 m _ 2 + 1 1 a,* ^ (0V j ) ii 



m-2,1 j 



Therefore, 



lirxi-E 



lim 

ej.0 



1 

2s~ 



d x id x 



^ip j )(y + x)y k (27Ts)- d/2 e- 1 ^ dyds 



- j d x ,d xJ (^^)(y + x)y k (27rs) 



(by (|3.14p and Fubini Theorem) 



^){x + y)y k 



(2TT S y d/2 e- 1 ^r dsdy 



dyds 



(3.15) 



2s 



= C d / d x id x j{(t> i tl) : >)(x + y)- ns dy 

jR d \u\ 
= d xk (-A)- 1 d x ,d x3 (^)(x) 

which coincides with the convolution representation of the operator V(— A) -1 described in [30l Page 31]. 
Hence, BSDE (|3.8|) belongs to class £ and by (|3. 14[) . Y £ C(R d ) on account of the continuity of the 
translation operator on L p (M. d ), p 6 [1, oo). Moreover, we have 



£(x) = V(-A)~ 1 ^ 1 ^ J (^(x)^'(a;)) = IkaE 



Vx e 



which completes the proof. 



□ 



Remark 3.3. We have Yq g C(R d ) from (|3. 14p and the continuity of the translation operator on L p (M. d ) 
for p 6 [1, oo). However, in this work we would rather use more properties of Yq basing on harmonic 
analysis. In fact, the operator P := I — VA _1 div is the well-known Leray- Hodge projection, where I is the 
identity operator. Note that the singular integral operator P (see [3111 GUI) is a bounded transformation 
in H nq for q G (1, oo) and neZ. Define P^ := I P. We have in Lemma l3~2l that £ = -P ± (div(0(g) ip)). 

Remark 3.4. For the defined £, there exists a scalar- valued fucntion rj E H m (M. d ;M.) such that £ = 
V77. Indeed, we may take r]{x) =: (— A)~ 1 d x id x j ((^>y)(;r) which, by the second order Elliptic partial 
differential equation theory (see [19]), lies in H m (R d ; R). 

Define the following operator 9: for any map g defined on [0,oo) x [0,oo) x R d , 

6 g (t, x) := git, i, x), (t, x) £ [0, 00) x R d . 



Our main result is stated as follows. 
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Theorem 3.3. Let v > 0,G G H™, and f G L 2 (0, T; i?™" 1 ) m > d/2. Tften our 

(|1.3p admits one and only one local H m -solution (X,Y, Z,Yq) on some time interval (Tq,T] with 9y G 
Q OC ((T , T]; H™) n L 2 oc ( T o, T; ff™ +1 ) and 9 Z G C ;oc ((T , T]; ff™- 1 ) n ^L( T o, T; ff™), uAere T depends 
on ^, T, ||G|j m and \\f\\L 2 (o,T;H m - 1 )- Moreover, there hold the following representations 

e z (t, •) = V9 Y (t, •), Y s (t, •) = 9 Y (s, X s (t, •)) and Z.(t, •) := 9 z (s, X,(t, •)), T < t < s < T, 

and there exists some scalar-valued function p such that Vp = Y and (By, 9 Zl p) satisfies 

d Y {r,X r (t,x))=G(X T {t,x))+ f [f(s,X s {t,x))+Vp{s,X s {t,x))]ds 

J r 

T 

9 z {s,X s (t,x))dW s , T < t < r < T, a.e.x G R d , a.s.. 

In addition, (9y ,p) coincides with the unique strong solution to Navier-Stokes equation: 
d t 9y + %A9y + {9y -V)9y + Vp + / = 0, T < t < T; 

V-9 Y = 0, 9 Y {T) = G. {Ab) 

Remark 3.5. In Theorem 13.31 we only have the connections for the strong solutions. For the case of 
to < d/2, we could not show that 9y(t, •) takes values in W^^° , and thus we do not know whether the 
forward SDE of FBSDS (|1.3[) is well-posed or whether the assertions of Lemma \3. II are still true. 

4 Existence and uniqueness 

4.1 Auxiliary results 

Consider the following coupled FBSDE: 

dX s {t,x) = [b(s,X s (t,x)) +aY s (t,x)] ds + ^dW s , s G [t,T]; 
X t (t,x)=x; 

-dY s (t,x)=</>(s,X s (t,x))ds- v ^Z s (t,x)dW s , s£[t,T]; 
Y T (t,x)=iP(X T (t,x)), 

where v > and a are constants. 

We define a bounded and a locally bounded solutions to FBSDE (|4.ip . 



Definition 4.1. We say (X,Y,Z) is a bounded solution (locally bounded solution, respectively) of 
FBSDE (gU) if for each t G [0,T] (t G (T ,T\ for some T G (0,T), respectively) and almost every 
X G M d , 

(JT.(t,aj),y:(t,aj),Z.(i,a:)) € 5 2 ([t,T];M d ) x S 2 ([t,T];M d ) x L^(t, T; M dxrf ) 

such that the forward SDE and BSDE on [0,T] (any [Ti,T] C (T ,T), respectively) hold almost surely. 

Lemma 4.1. Let b,c/> G C^°(E d+1 ) and V G C^°(R rf ). T/ien /or eac/i integer to > d/2, FBSDE (|4TTj) 
admits a unique locally bounded solution (X,Y,Z) on some time interval (r, T] wii/i T — r (V G [0,T)J 
continuously depending on \\4>\\L 2 (o.T-H m - 1 )i \\b\\c([o,T];H m )i HV'llm, v and a. Moreover, we have 

9y G C ioc ((r,T] ;j ff m ) n L 2 1oc (t,T; H m+1 ), (4.2) 

and /or eac/i t G (t, T] , almost all x G R d and all r G [t, T] 

9 Y (r,X r {t,x))=^{X T [t,x))+ f ${s,X s {t,x))ds-^D { 9 z {s,X s (t,x)) dW s , a.s. (4.3) 



Z (M) = Wy (*.*)> Y r (i,:c) =6»y(r,X r (t,x)), ^(t,a:) = ^(r,X r (t,a!)), a.s.. (4.4) 
In particular, if a — 0, we are allowed to take Tq = and we /iaue 

0y G C([0, T]; iJ m ) n L 2 (0, T; H m+1 ) n C*°°([0, T] x R d ). 
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Lemma |4. II might exist elsewhere, but we have not found it. For the reader's convenience, a proof is 
sketched in the appendix. 

Remark 4.1. In Lemma |4~TI for each T <E (r, T) and each T" € (T ,T], consider the following FBSDS: 

dX s {T 0l x) = [b(s,X a (T ,x)) + a9 Y (s,X s (T ,x))} ds + ^dW s , s € [T ,T']; 
X To (T ,x) = x; 

-dY s (T ,x) = <l>(s,X s {To,x))ds - V^Z s (T ,x)dW s , s£ [T ,T f ]; 
{ Y t ,(To,x)=0y(T',Xt'(T o ,x)). 

It has the same unique locally bounded solution (X.(Tq, x), Y. (To, x), Z.(Tq, x)) as FBSDS (|4.1|) . Define 
the following new equivalent probability measure: 



In view of Girsanov Theorem, the process 



expf-^ I {b + aOY^XsiT^x^dWs^^v- 1 I ^ \{b + a6 Y )(s, X S (T , x))\ 2 ds)dF. 



W r :=v—*\ {b + a6 Y ){s,X s (T ,x))ds + W r , re[T ,T] 
Jt 

is a standard Brownial motion under Q T °' X and we have 
f dX s (T ,x) = V^dW s , s e [T ,T']; 
X To (T Q ,x) = x; 

-dY s (T Q ,x)= [cj ) { S ,X s (T ,x))+e z (b + a9 Y )( S ,X s (To,x))]d S -V^Z s {To,x)dW s , s € [T ,T'}; 
{ Y T ,(T ,x) = u(T',X T ,(T ,x)), 

where we have used the expression for Z in (|4.4j) . For the unique solution (X'.(Tq, x), Y.'(Tq, x)), Z[(Tq, x) 
to the following FBSDS: 

f dX' s {T ,x) = V^dW s , s E [T ,T'}; 
x t„( t o,x) = x; 

-dY^T ,x) = I m [4> + 9 z (b + ae Y )](s,X' s {T ,x))ds-V^Z' s (T ,x)dW s , sG [T ,T'] ; 
Y^,(T ,x) = I m u(T\X T >(T ,x)), 

it is not difficult for us to show that for almost every x € K d and all s£ [T , T'] 

Y^(T Q ,x) = I m e Y (s,X' s (T ,x)) and Z' s (T Q ,x) = I m 9 z (s, X' S (T , x)), a.s., 

and from Ito's formula, we have 



r 



\Y± Q {T 0l x)\ 2 +v \Z' s (T ,x)\ 2 ds 



T' 



E QTo ,,[< n'(T ,x), I m <t>(s,X' s (T ,x)) + I m {d z (b + a6 Y )(s,X' s (T ,x))} > ]ds 

+ E 9 t „ [\I m 9 Y (T',X^(T ,x))\ 2 ] , a.e.x e R d . 
Finally, integrating with respect to x on both sides of the last equality, we obtain the energy equality: 



\\9 Y (T )f m + v 



T' 



h{s)\\ 2 m ds 



T 



Note that 



\\9 Y (T')\\ 2 m + 2 / (4>(s) + 9 z (b + a9 Y )(s),9 Y (s)) m ^, m+1 ds. 

J T 



IIM*)llm = P Z (t)\\ 2 m-l + IM*)llo < Pz(t)fm-1 + WMl-V 
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Remark 4.2. From Lemma [4.11 we see 9 Y and 9z are all deterministic functions. Therefore, for each 
semimartingale {X' s (t,x), s £ [t,T]} of the form 

X' s (t, x) = x + J <p r (t, x)dr + J^ s/vdW r , T <t<s<T 

with {ip s {t, x), s £ [t,T]} being bounded and predictable, it is interesting to understand (6 Y ,9z)(s, X' s (t, x)) 
in the FBSDE framework. Indeed, Analogous to the preceding remark, define the following equivalent 
probability measure: 



1 

2^ 



:= cxp J [(b + a9 Y ){s, X' s (t, x)) - <p.(t, x)] dW s 

\{b + a8 Y ){s,X' s {t,x)) - tp s (t,x)\ 2 ds 
Then in view of Girsanov theorem, there is a standard Brownian motion (W',<Q' X ) such that 

X' s (t, x) = x + J (b + a6 Y )(r, X' r (t, x)) dr + J JDdW' r , T <t<s<T. 
From Lemma |4.1[ we have 

6 Y (r, X' T {t, x)) = ^{X' T {t, x)) + J* X' s (t, x)) ds-V^f e z (s, X' s (t, x)) dW' s 

= ^{X' T {t,x)) + J {[<t>(s,X' s (t,x)) + 6 z (b + a9 Y )](s,X' s (t,x)) - Z<p a (t,x)} ds 
9 z (s,X' s (t,x))dW s , t<r<T. 



Lemma 4.2. Let b £ C([T Q ,T}; H m ) n L 2 {T ,T;H m+1 ),(t> £ L 2 (T a ,T; H" 1 ^ 1 ), and tp £ H m , with 
m > d/2. Then, for each t £ [Tq,T) and almost every x £ R d , the following FBSDE: 

dX s (t,x) = b(s,X s (t,x))ds + y/vdW s , T — e <t < s <T; 
Xt(t, x) — x\ 

r- (4-5) 

-dY a (t,x) = 0(s,X s (t,x))ds- V^Z s (t,x)dW s , s£ [t,T]; 
Y T (t,x)=iP(X T (t,x)) 

admits one and only one solution 

(X.(t,x),Y.(t,x),Z.(t,x)) £ S 2 (T ,T;R d ) x S 2 {T ,T;R d ) x L%{T ,T;R d ), 

and for this solution {X,Y,Z), there hold (|4.2p , (|4.3[) and (|4.4[) with r being replaced with T therein. 

Proof. Form > d/2, H m ^ C°' S , H m+1 <-+ C 1 ' 6 . By Theorems 3.4.1 and 4.5.1 of .27 , the forward SDE is 
well posed for each (t, x) £ [To, T] x R d and defines a stochastic flow of homeomorphisms. Moreover, from 
Lemma \3. II and Remark 13.11 the backward SDE is also well posed for every x £ R d /F t with Lebesgue's 
measure of F t being zero. Therefore, for each (t,x) £ [T ,T] x (R d /F t ), FBSDE (|4.5jl has unique solution 

(X.(t,x),Y.(t,x),Z.(t,x)) £ S 2 (T ,T;R d ) x S 2 {T Q ,T;R d ) x L%{T ,T;R d ). 

For each (t,x) £ [T ,T) x (R d /F t ), define the following equivalent probability measure: 

dQ Lx :=exp (^-^J f b{s,X s (t,x))dW s - ^ jf \b(s, X s (t, x))\ 2 ds^j dP. 
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Then there is a standard brownian motion (W , Q ttX ) such that FBSDE (|4.5|) is written into the following 
form: 



(4.6) 



' dX s (t, x) = VvdWs, T - e <t < s < T; 
X t (t,x) = x\ 

-dY s {t,x) = 4>(s,X s (t,x)) + Z s (t,x)b(s,X s (t,x))ds - ^Z s (t, x) dW' s , s € [*, T] ; 
Y T (t,x)=xP(X T (t,x)). 

Choose a sequence ,f Z 1 } C C" c 00 (M rf+1 ) x C c °°(R rf+1 ) x C~(R d ) satisfying 

lini^ - b\\ C ([t,T];H^) + \\ bn - b\\ L 2([t,T];H™+i) + U n - <£|U=(t,T;iT"-i) + \W ~ V'llff'") = 0. 

Let (X,Y n ,Z n ) be the unique solution of FBSDE gU) with (b,(/),ijj) being replaced with (b n , (j> n , ij) n ). 
Then for each n, we have 

<9yn g C{[t,T]-H m ) nL 2 (i,T;iJ m+1 ), 
and for each a; € R d and all i < r < T, 

9 Y n (r, X r (t, a:)) = V" (Xr(i, x))+ [ (<f> n + Zn b n ) (s, X s (t, x))ds-^[ 0z~ {a, X s (t, x)) dW' s , 



0z»(t,a:) =Vd Y n(t,x),Y?{t,x)=9 Y n{r,X r (t,x)),Z?(t,x) = 6zn(r,X r (t,x)). 
On the one hand, using Ito's formula, we have, for each (s, x) € [t, T] x M. d /Ft 

/■T 

|r s n (i,a;) -r s (i,a;)| 2 + ^ y Iz; 1 ^,^) - ^r(i,x)| 2 dr 

= |^(X T (t, x)) - ip n (X T (t, x))\ 2 ~ J 2u(Y r n (t, x) - Y r (t, x), Z?(t, x) - Z r (t, x)) dW' r 

+ J 2(Y; l (t,x) - Y r (t,x), (<f> n - <P)(r,X r {t,x)) + Z?(t,x)b n (r,X r (t,x)) 
- Z r (t, x)b(r, X r (t, x))) dr, 
and therefore, using the BDG inequality, 



Ent.t 



< Em,* 



1 



sup \Y?(t,x) -Y T (t,x)\ 2 + is I \Z?(t,x)-Z r (t,x)\'dr 

re[s,T] 

mx T (t, X ))-r(x T (t,x))\ 2 



+ - sup \Y?(t,x)-Y T (t,x)\ J + C / \Z?(t,x)-Z r (t,x)\ dr 

^ r£[s,T] Js 

+ C f \Y r n (t,x) - Y r {t,x)\{\<j> n - <b\(r,X r (t,x)) + \\b n - b\\ c(%T] . Hm) \Z r (t,x)\ 

J S 

+ \Z?(t,x)- Z r (t,x)\)dr 
with the constants Cs being independent of n. Combining (|4.7p and (|4.8p . we have 



Em, 



< CEnt,x 



sup \Y?(t,x)-Y T {t,x)\ 2 + v / |Z"(i,:z) -Z r (t,x)\ dr 

re[t,T] Jt 
rT 

(\cf> n (r,X r (t,x)) - <Xr,X r (i,x))| 2 + \\b n - b\\ 2 c([t<T] . Hm) \Z r {t,x)\ 2 ) dr 



(4.7) 



(4.8) 



+ C\\il) - ip n \\ 2 Hm — >0, asn->oo, xeR d /F t . 



(4.9) 
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On the other hand, setting 

(Y nk , Z nk ,b nk ,cb nk , i> nk ):= (Y n - Y k , Z n - Z k , b n - b k , 0" - <f> k , ^ - ^ k ). 
For n, k g Z + , we have from Remarks I2.1l and l4.1l that 



\\6 Ynk (s)r m + v J \\9 Znk (r)\\t n dr 

= ll^nfcllm + / 2(4>nk(r) + Z nb n k{r) + 9 Znk b k (r), 9 Yn h(r)) m ^ 1>m+1 dr 



<Unkf m + ^ J (\\9 Z nH(r)\\ 2 m + \\e Y n k (r)\\ 2 m )dr 

+ C(u) ( f T \\b nk {r)\\l\\8Y»(r)\\ldr+ f \Mr)\\ 2 J\9 Ynk (r)\\ 2 m dr + / |K,. (r)||f„_ L dr 



<C{U nk \\ 2 m + - J {\\e znk (r)r m + \\9 Ynk (r)rjdr 

T pT 

2 i II J, . /„M|2 \ j_ i / \\q . /„\||2 



(IIM»-)ll^ + ll^fc(r)|lS,-i) ir + J \\9Y^{r)t m dr' 
where we have used the the following priori estimate by taking (b k , cf> k , ip k ) = in the above, 

\\8y<>\\c([To,T]:H™) + \\8Y™\\L 2 (T ,T;H m + 1 ) - @, 

with the constant C being independent of n. Thus, 



sup \\9yn>(8)C + u f 1 \\9 zn *(r)\\ 2 m dr 

se[T ,T] J T 



(4.10) 



< C \\\^nkf m + J t (IIMOIIm + \\4>nk(r)\\ 2 m -i) drj -^0asn,k^^. 
Combining (|4.9j) and (|4.10p . we have 

lim (\\9 Y k - 6y\\c{\T ,T];H m ) + \\0z k - 6z\\l%(T ,T;H™)) = 0, 

with 9z = V0y, and furthermore, by taking limits, we can deduce that (|4.2|) . (|4.3[) and (|4.4j) hold. We 
complete the proof. □ 

Proposition 4.3. Assume that ip e H m , b G C([0, T]\H m ) n £ 2 (0, T\ H' m+1 ), and <j) G L 2 (0, T; if™ -1 ) 
with m > d/2. Then FBSDE (|4.1I) admits a unique locally bounded solution in some time interval (Tq,T] 
with To depending on \\ip\\m> \\b\\c([o,T];H m )> \\<t l \\L 2 (o.T;H m - 1 ), a , v an d T. If a = 0, To = 0. Moreover, 

9y G C loc ((T ,T};H m )nL 2 oc (To,T;H m+1 ) and 9 Z e C loc ((T ,T}; H m ~ l ) n L 2 OC (T ,T; H m ) 

satisfy ()4.3[) and for any t S (To,T], i/iere /ioW the following energy equality: 

pY{t)\\ 2 m + v j \\6 z {s)\\ 2 m d S 

= U\\ 2 m + 2 (0 z (b + a9 Y )(s),0Y(s)) m -i, m +ids + 2 J (<£(s), 6» r (s)> ro _i, m+1 ds 
and for almost every x £ M. d and all s € [t,T], 



(4.11) 



(t,x) = V9 Y (t,x),Y s (t,x) = 9 Y {s,X s (t,x)), Z s (t,x) = 9 z (s,X s (t,x)), a.s.. (4.12) 



In addition, if m > + 1, t/te locally bounded solution on the time interval (Tq,T] is the unique 
local solution as well. 
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Proof. By Lemma 14.21 it remains for us to consider the case a/0. 

Step 1. We shall prove the existence of the solution. Choose a sequence (b n , </>„, tp n ) £ C£°(R 1+ ) x 
C*o°(K 1+d ) x c™(R d ) guch that 

Jirn ) {||6„ - b\\ C ([0,T];H™) + \\ b n - &||l,2(0,T;fr™+ 1 ) + \\<Pn ~ <P\ | L 2 (0,T;H™~ 1 ) + llV'n ~ ^|U} = 0, 
\\bn\\c([0,T];H™) < C||^llc([0,T];_ff m ) 5 1 1 </>n 1 1 L 2 (O.T-H^- 1 ) < C|| 011 L 2 (0,T;ff™- 1 ) , llVVillm < C'UV'llm, 

and ||&n||L2(o,T;/f m + 1 ) < C||^llL 2 (o,T;ff m + 1 )i where C is a universal constant being independent of n. By 
Lemma B~T1 for each n, FBSDE (|4. 1[) with (6, 0, ?/;) being replaced with (6 n , </>„, admits a unique local 
solution (X n ,Y n , Z n ) on some time interval (r, T] such that (9 Y -n.,0z") satisfies (|4.3[) associated with 
(b n ,<t> n ,ip n ). As T -t continuously depends on ||0 n |U 2 (o,T ; H"'-i) 5 ||Vv»||to, \\K\\c([o,T];H™), v and a, we 
can choose a uniform t for all n £ Z. Moreover, we have 

9 Z n(t,x)=V9 Y n(t,x),Y s n (t,x) = 6 Y n{s,X2{t,x)) and Z?(t,x) = 0z»(a,X?(t,x)), 

and by Remark 1 4. 11 

\\0Y4s)f m + v f Pz~{r)\\ 2 m dr 

J S 

= \\i>n\\m + J 2(9 Z n(b n + a9 Y n)(r),e Y «(r)) m - ltm +idr + J 2{(</> n (r), 6 Y n(r)) m _ hm+1 dr 



< CM| £(\\ b (r)\\ 2 L ^ , T , Hm) + l) \\e Y n{r)t m dr + o? Jj6 Y n (r)C dr 

2 J„ I / ||/)_ ^„M|2 ; , /"f||„/.||2 



Therefore, 



II^WIIm-idrK 2 / ll«z-(r)||^dr + C|^| 



T /.T 

2 j„ , „2 / Mfl / mi4 



< C{v,<j>,ijj)+C(v,b) I y ||^(r)||^ dr + c^y ||0y»(r)||^dr 

< C(0,V,^&,«,T) + C(^6)a 2 / \\e Yn (r)f m dr. 



In conclusion, there is a constant To = To(</>, ^, 6, a, T) £ [r, T) such that for Co := C{4>, ip, v 1 b, a, T), 
sup ||^n(r)||^ + i./ T ||^n(r)||^dr<- ^° -, V S e(T ,T]. (4.13) 



Set 

(X nk , y ™ fe , Z nk , 6„ fc , cf> nk , V„ fc ) : = (X" - X fc , Y " - y fc , Z*» - Z fe , 6 n - 6 fc , 0„ - fe , - Vfc ) ■ 



1G 



Then for each hxed e G (0, T - T ), we have, for any s E (T Q + e, T) 



\\9 Ynk (s)\\ 2 m + v ( T \\9 Znk (r)\\ 2 m dr 



ll^nfc|| m + J 2{8z*>b nk (r) + 6 Z nkbk(r), 9 Y nk(r)) m ^x,m+idr 

+ J 2a(9 Z nO Y nk(r) + 9 z „ k 9 Y k(r), 9 Y nk(r)) m - ltTn+1 dr + J 2(<f> nk (r), 9 Y nk(r)) m -i, m+1 dr 



< Unk\t + ^ I {\\9 Znh {T)\\ 2 m + \\9 Yn u{r)\\ 2 Jdr 



+ C(«/) | jT \\hnk(r)\\ 2 m \\9 Y n(r)\\ 2 m dr + jT \\b k (r)f m \\9 Y n k (r)f m dr 

< W^nkWl + (P znh (r)\\ 2 m + \\9 Ynk (r)\\ 2 Jdr 

+ C[ f (IIMOIIm + dr+ j T \\9 Y ,Ar)\\ 2 m dr\ . 



Consequently, 



sup \\9 Ynk (s)\\ 2 m + v [ \\9 Zn u{r)f m dr 

se[T +e,T] JT +s 

< C (Unkfm + £ (\\b nk (r)f m + ||<^ fe (r)ll™-i) dr J — > as n, k 



(4.14) 



where the constant C is independent of n and k. Define 9 Z n '■= V0y> for n G Z + . Then {(9 Y ™, 9 Z ™), n € 
Z+} is a Cauchy sequence in C([T + e, T]\H m ) x L 2 {T + e, T; H m ), whose limit is denoted by (£ VQ. 
By taking limits and in view of the arbitrariness of e, we check that the pair (£, V£) satisfies the energy 
equality (|4.11l) on (T ,T] with (9 Y ,9 Z ) being replaced with (£,V£). Moreover, through a bootstrap 
argument, we can extend the existing interval to a maximal one still denoted by (To, T] with To depending 
on ||V>|L, \\b\\c(lo,T];H^), \\4>\\L^(o,T;H—^),a,v and T. 
On the other hand, consider the following FBSDE: 

dX s (t,x) = [b(s,X s (t,x)) + aC(s,X s (t,x))] ds + ^/udW s , se[t,T\; 
XAt, x) = x; 

V (4.15) 

-dY s (t,x) = (f>(s,X s (t,x))ds- y/uZ s (t,x)dW s , sG [t,T]; 

Y T (t,x) = 1 p(X T (t,x)), 

which admits a unique local solution (X, Y, Z) on (T , T]. From Lemma[01 we have (|4~21) . (|4"3)) and (|4~4")l . 
Letting fc ->• oo in P~T4"1) . we have Cfos) = 0y(*,a:) and VC(i, a;) = 9 z {t,x) for a.e. (t, a:) G (T ,T] x M d . 
Furthermore, from Lemma FOl we deduce that the triple (X s (t,x)X(s 7 X s (t 7 x)),V((s,X s (t,x))) solves 
FBSDE (|4.1[) and satisfies all the assertions of this proposition except the uniqueness, which is left to the 
next step. 

Step 2. We now verify the uniqueness. Let (X, Y, Z) be any locally bounded solution of (|4.1j) on 
(Tq, T]. For each t G (To, T] and almost every x G R d , define the following equivalent probability measure: 



:=cxp (^-^J t [b(s,X s (t,x)) + aY s (t,x)} dW s - ^ |6(s, X a (t, z)) + a y s (t, a;)| 2 rfs^J 



dl 
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Then FBSDE (03} reads 

dX s {t,x) = ^dW' s , s G [t,T]; 
X t (t,x) = x; 

-dY s {t,x)= (f>(s,X s (t,x)) + Z s {t,x)(b(s,X s (t,x)) + aY s (t,x)) ds-^Z s {t,x)dW^, s€ [t,T 

Y T (t,x)=^(X T (t,x)), 

where (W',Q t:X ) is a standard Brownian motion. 
Define 

Y?(t, ■) = 0y„( S ,X s (t, •)) and Z?(t, •) = 6 Z n(s,X s (t, ■)). 
As m > d/2 and iJ m ^ C M (M d ) for some 5 G (0, 1), there is a constant N l such that 



sup sup 

n Vse[t,T],x£ 



(4.16) 



y»(s,a;)| + / sup Z "(s,x) < AT*. 



By Remark |4. 2 [ we have for almost all x &M. d , 

Y s n (t,x) =il> n {X T (t t x))-W ( Z?(t,x)dW; 



+ / [Z?(t,x)(b n (r,X r (t,x)) +aY r n (t,x)) + (j) n {r, X r (t, x))] dr, t<s<T. 

J s 

Then Ito's formula yields 

\Y s n (t,x)-Y s (t,x)\ 2 + v [ \Z?(t,x)- Z r (t,x)\ 2 dr 



\tp n (X T (t,x)) - ^j(X T (t,x))\ 2 + 2 J (Y r n (t,x)-Y r (t,x), 

Z?(t,x)(b n (r,X r (t,x))+aY?(t,x)) - Z r (t,x)(b(r,X r (t,x)) + aY r (t,x)) 

+ <t> n {r,X r {t,x)) ~ cfi{r,X r (t,x)))dr-2 [ (Y r n {t,x) - Y r (t,x), (Z? - Z r )(t,x) dW r '). 



Note that both (Y n (i, x), Z n (t, x)) and {{Y.{t,x),Z.{t,x))) belong to S 2 ([t,T];R d ) x L%(t,T;R d x R d ) 
and moreover, there exists a constant K*' x such that sup sg r t T i \Y s (t, x)\ < K f ' x , a.s.. Then, we have 



Ent.x 



\Y s n {t,x)-Y s (t,x)\ 2 + v / \Z?(t,x)-Z r (t,x)\ 2 dr 



(Y r n (t,x)-Y r (t,x), 



= Eqt.x Un(X T {t, X)) - 4>(X T (t, X))\ 2 } + 2E®t,: 

Z r r l (i, x) (b n (r, X r {t, x) ) + aY r n (i, x)) - Z r (t, x) (b(r, X r (f , x)) + aY r [t,x))+(/> n (r, X r (t,x)) 
-<t>(r,X r (t,x)))dr 



v 



\Z 7 r l {t 1 x)-Z r (t,x)\ 2 dr\ +C{m,d)U n -ip\\; 

l-T 



+ C{v, a)Eqt 



(\</> n (r,X r (t,x)) - ^{r,X r (t,x))\ 2 + N^bn - b\\ 2 c([0 T] . Hm) 
\Y; l (t,x) -Y r (t,x)\ 2 (l + \\b(r)\\ 2 n + \\e Zn (r)\\ c{Rd) + \Y r (t, x)\ 2 )) dr 



<\E q 



CEnt,* 



\Z?(t 1 x)-Z r (t,x)\ 2 dr\ +C{m,d)U n -nL 

T 

{\<j) n {r,X r [t,x)) - cj>(r,X r (t,x))\ 2 + \\b n - 6||c([o,t] ; ^ 
+ \Y r n (t, x) - Y r (t,x)\ 2 {l + \\6 z ™(r)\\c { R«))) dr 
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Then, from Gronwall Inequality, we have 



sup E Q 

se[t,T] 



\Y s n (t,x)-Y s (t,x)f 



< C 



C{[0,T];H" 



En 



\Z™{t,x)-Z s (t,x)\ 2 ds 



(\4> n {r,X r {t,x)) - ^{r,X r {t,x))\ 2 ) dr] } 



0, 



where the constant C depends only on iV', K t,x , T, ||&||c([o,T];.H' m ): T, m, d, v and a, and is independent of 
n. 

Thus, in view of (|4.16p . we conclude that 

Y s (t, x) = C(s, X a (t, x)) and Z B (t, x) = VC(s, X a (i, x)). 

Therefore, any locally bounded solution of FBSDE ()4.1j) on (Xb,T] must have the form described as 
above. Now, let (X, Y, Z) and (X, Y, Z) be any two locally bounded solutions of FBSDE (ET~Tj) on (T , T]. 
By previous argument we have 



Y s (t,x) = C(s,X s (t,x)), Z s (t,x) = V{(a,X,(t,x)), 
Y s (t,x) = ((s,X s (t,x)), Z s (t,x) = V((s,X s (t,x)). 



(4.17) 



Hence X. (t, x) and X. (t, x) satisfy the same forward SDE with the same initial value. Thus we must have 
X = X, a.s., which in turn shows that (Y, Z) = (Y, Z), a.s.. 

Step 3. For m > d/2 + 1, to prove that the unique locally bounded solution on (T ,T] is also the 
unique local solution, it is sufficient to prove that the locally bounded solution constructed in Step 1 



(X s (t, x), C(s, X s (t, x)), VC(s, X s (t, x))) To<t 



<s<T 



is the unique local solution to FBSDE (|4.1[) on (T ,T] as well. 

Since m > d/2 + l and H m ~ l ^ C°' S (R d ) for some 5 € (0, 1), our BSDE is well-posed for each x G M d . 
Let (X, Y, Z) be any solution of flUT) on [i, T] with 4 G (T , T]. For every a; G K d , define 



r; i (t,x) = 0y»(s,X s (t,a:)) and = 0z»(s, a;)). 

By Remark |4. 2 [ we have 

Y s n {t, X) = 1> n (X T (t, X)) - y/V I Z?(t, X) dW r 



(4.18) 



[Zl!(t,x)(b n (r,X r (t,x))+aY 7 n (t,x) - b(r,X r (t,x)) - aY r (t,x)) 
+ 4> n (r,X r (t,x))] dr, t <s <T, a.e. x G R d . 



Then Ito's formula yields 



\Y?(t,x)-Y.(t,x)\ J + u / \Z?(t,x)-Z r (t,x)\'dr 



\<P n (X T (t,x)) - ^(X T (t,x))\ 2 + 2 J (Y r n (t,x)~Y r (t,x), 

Z?(t,x)(b n {r,X r {t,x)) + aY r n {t,x) - b(r,X r {t,x)) - aY r (t,x)) 

+ <l> n (r > X r (t,x))-cfi(r,X r (t,x)))dr-2 [ (Y r n (t, x) - Y r (t, x), (Z? - Z r )(t,x) dW r 



Note that both (Y n (t, x), Z n (t, x)) and {{Y.(t,x),Z.{t,x))) belong to S 2 ([t, T];R d ) x L%(t,T;R d x R d ) 
and moreover that there exists a positive constant K t such that 



sup sup \Z?(t,x)\ < C\\e Z n\\ ci[t}T] . Hm - 1) < K t , a.s., Vx G R d . 

n se[t,T] 
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Thus, we have 



\Y s n (t,x)-Y s (t,x)\ 2 + v / \Z?(t,x)-Z r {t,x)\ 2 dr 



E 



E[\*p n (X T (t, x)) - 4>{X T {t,x))\ 2 ] + 2E 



(Y r n (t,x)~Y r (t,x), 



Z?(t, x) (b n (r, X r (t,x)) + aY, n {t,x)~ b(r, X r (t , x) ) - aY r (t,x)) + <j> n (r, X r (t,x)) 
— (j)(r,X r (t,x)))dr 

< C(K t ,T, a, v, d, m)E J \Y r n (t, x) - Y r (t, x)\(\Y r n (t, x) - Y r (t, x)\ + \\b(r) - b n (r)\\, 
+ \\cf> n {r) - 4>{r)\\ m -i) dr + C\\i/> n - i>f m . 
Using Gronwall inequality, we have 

sup E\\Y s n (t,x)-Y s (t,x)\ 2 
se[t,T] L 



vE 



\Z^(t,x)~Z s (t,x)\ 2 ds 
"iff™- 1 )} > 



(4.19) 



2 



0. 



< C{\\ip n - ipW^ + \\b - b n \\ 2 c[[G T] . Hm) + \\(j) n 
Thus, in view of f|4. 18[) . we conclude that 

Y s {t,x) = ((s,X s (t,x)) and Z s {t,x) = VC(s, X s (t, x)). 
Analogous to the arguments of Step 2, we verify that the triple constructed in Step 1 

(X s (t, x),C(s, X s (t, x)), VC(s, X s (t, x))) To<t < s < T 
is the unique local solution to FBSDE (I4.1j) on (To, T] as well. The proof is complete. 



□ 



Remark 4.3. In Step 2, the equation (|4.3I) plays a crucial role in the proof of uniqueness. Indeed, 
FBSDE (|4.1[) is usually associated to the deterministic vector field By which satisfies (|4.3[) instead of 
being a classical solution to some parabolic PDE in [29]. Note that equation (|4.3[) is probabilistic and 
that according to Lemma 13. 1[ it makes sense for our Sobolev coefficients. Therefore, our method here 
helps to probabilistically solve more general coupled FBSDEs. 

On the other hand, in view of the whole proof of Proposition ^. 31 we have 



T < V 



T — 



1 



a 2 C C(v, b) 



and also, if Ta 2 CoC(v, b) < 1 in (|4.13[) . then the local bounded solution is actually a global solution on 
the whole interval [0,T]. 

Remark 4.4. If F (-, •) of FBSDS JO]) lies in L 2 (T , T; H m ~ l ), then 9 Y of (fOJ) is deterministic and 
belongs to L 2 (Tq,T; H m+1 ). Therefore, by Lemma [3.11 and Proposition 14.31 Definition 13.21 and Remark 
make senses. 



From Proposition I4.3[ we have the following characterization of an _ff m -solution to FBSDS (|1.3|l for 
to > d/2, whose proof is omitted. 

Corollary 4.4. Under assumptions of Theorem \3.3l (X, Y,Z, Yq) is an (local, respectively) H m -solution 
of FBSDS (|1.3[) (on some time interval (To,T] ; respectively) if and only if (X, Y,Z, Yq) is a solution of 
FBSDS CL3]) with By e C{[0,T]:H m ) (By € C loc {{T 0l T}; H m ), respectively). 
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4.2 Proof of Theorem EH 



First, for each v € C{[Q,T];Hf) n L 2 (0,T; H m+1 ),( G C([0,T];JJ m ) n L 2 (0,T; H m+1 ), consider the 
following FBSDS: 



-dF s (*,x 
-dY s (t,x 



v(s,X s (t,x))ds + y/vdW a , s G [t,T] 



/(«,Jf,(*,as))+F (*,^«(*,a!))J da- ^uZ a (t,x)dW a , ae[t,T\; 
G(X T (t,x)); 

^C(t, x + b s ) (sl - si) - m ) B f 

- dM s , s G (0,oo); 



(4.20) 



k r 00 (*,s) = o. 

By BSDE theory (see [25l|33]) and Lemma[3J2 FBSDS (lL20| admits a unique solution (X v ' c , Y v ^, Z V '^Y^ C ) 
and in view of Lemma 13.21 and Remark 13.31 there holds that 



Y (t,x) = -P X (div(v(t,x) <g> ((t,x))) = -P X ((v(t,x) ■ V)C(i,x)), 

where we have used the fact that div(v) = 0. From Proposition 14.31 and Remark 14.31 we conclude that 
(0y,.c,0z.,c) € C([0,T];H™) x L 2 (0,T;H m ). 

For any & G C([0, T]; _ff m ) n L 2 (Q, T; i/ m+1 ), i = 1,2, set 



(W y „,c,^,c,<5C) 
Then by Proposition ^. 3[ we have 



, C1-C2). 



\\66 Y .. ( (8)\\ 2 m + v J \\86 z .. ( {r)\\ 2 m dr 

2(59 z *, c v(r), %„{ (r)) TO _i, m+ i dr - / 2(P ± ((i- • V)5(){r), 86 Y *.i (r)) m -i, m +i dr 

J S 

< ^j\\\S0 Y ^(r)\\l + \\5e z ,. x (r)\\ 2 m )dr 

T fT \ 

2 WXO _/"^M|2 j„ 1 / ll.,/'_M|2 ||rA('„M|2 



+ |Kr)0*Wr)||jUr+ / II <->ll7, II'M <->ll7, '/<'• 

Using Gronwall inequality, we obtain 



sup ||^..c(«)llm+ / \\S9 z ,x{rW m dr<C{T-t)\\5aoat,nB^) 
se[t,T] ' 



(4.21) 



with the constant C depending on v, ||v||c([o, T] ; £T m ) and T. Then, by the contraction mapping principle 
we can choose a small enough positive constant e < T depending only on ||v||c([o, T];£T m ) and T, such 
that there exists a unique function ( E C([T — e,T]; H m ) satisfying 

(<W,<W) = (C,VC) mC([T-e,T];H m ) x L 2 (T - e,T;H m ). 

Then by Lemmas 13.21 and 14.21 we have for almost all x G R d , 

6 zv ,z(t,x) = V6 Y ^(t,x), Y r ^(t,x) = 6 Y ^(r,X r (t,x)), Z v /(t,x) = z „, c -(r,* r (i, x)), a.s., 



9 YV ,z (r, X r (i, x)) =G(X T (t, x)) + / /(*, X,(t, a;)) - P x ((« • V) e YV ,- c )(s, X s (t, x)) 
>/ / ^„,<(s,X s (t,iE))dW s , a.s.. 



ds 
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For each (t,x) G [T — e,T) x R d , define the following equivalent probability measure: 



exp^--^jf v(s,X s (t,x))dW s - i- jf |6(s,X s (i,x))| 2 rfs^ dl 



Then FBSDS (j4^0|) reads 
rfX s (i, x 

-dr a (t,s 



= V^dWi, sG[t,T]; 



/ + 9 z v + Y 



(s,X s (t,x))ds- <JDZ s (t,x)dW' s , sG [t,T]; 



G(X T (t,x)); (4.22) 
^v%(t, x + B.) (BL - B|) (flrf - Bl) B f 
- dAf s , s G (0,oo); 

t r 00 (t,x) = o, 

where (W',Q t ' x ) is a standard Brownian motion. Then taking the divergence operator on both sides of 
the above third and fourth equalities, we conclude that 

div 6 YV - t (s,y) = 0, a.s. Vs G (t, T], a.e.y G M d . 

On the other hand, by Proposition 14. 31 and Lemma |2~21 we have 



\\Gf m + / 2<P((« • V)^,c)(r), V,sW>m-i,m+i dr + / 2(/(r), y ., f (r)) ro _i lW +i 



<||G||1 + C / |Kr)|U|^, f (r)||^dr+- / \\f(rW m -i dr 



+ ^ I (IIWOI& + IIWrC)* 



which together with the Gronwall inequality implies 

sup + £ l T \\e z ,Ar)\\ 2 m dr 



se[t,T] 



2 7 t 



< C(„,T) (||/||! 2(0 ,t;H— ) + IIGII^) e (ll»ll««*, a -i ! *» J +-)(=r-*>. 



(4.23) 



Hence, through a bootstrap argument, we conclude that there exists a unique function £ G C([0, T]; if™) 
satisfying (6> y „, f , 6 ZVti ) = (<f, V<f) in C([0, T]; B™) x L 2 {0, T; H™), and again by Lemma EH we conclude 
that 

(X V ,Y V ,Z V ,Y^) := (X V ' C ,Y V ^,Z V ^,Y^ X ) 
is the unique B m -solution of the following FBSDS: 

' dX s (t,x) = v(s,X s (t,x))ds + ^/vdW s , s G [t,T]; 
X t (t,x) = x; 

-dY s (t,x) = [f{s,X s (t,x)) +Y (s,X s (t,x)j\ ds-^Z s (t,x)dW s , sG [t,T]; 

Y T (t,x) = G{X T (t,x)); (4.24) 
27 



dY s (t, x) = —v l Y t J (t, a: + B s ) BL -51 B{- Bi B* S - J ds 



dM s , s G (0,oo); 



k y oc (i,x) = o. 
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Choose two positive real numbers R and s (e < T) whose values are to be determined later, and define 
Ur ■= { u e C{[T - e,T}; H™) D L 2 {T - s,T: H™ +1 ) : 

\H 2 C([T-e,T]- >HT ) + ^\\Vu\\ 2 L 2 {T _ e>T . H ™+i ) < R 2 }- 

For any v £ £/Jj, there holds the following estimate by (14. 23[) : 

sup \\e YV (s)\\ 2 m +^ f \\9 z *{s)\\ 2 m dr< C(i/ J ||/||i a(0)T . HW ._ 1) ,||G||S,,T)e Jte . 
Choosing R to be big enough and e to be small enough, we have 



(4.25) 



sup \\Q Y v{s)f m + V - f T \\6 z ,(r)\\ 2 n dr<R 2 . 

\T-e,T] z JT-e 



se[T-E,T] 

On the other hand, for any v\, t>2 £ f/fj, setting 

(SOy* , 69z* , Sv) := (0y»i - , 9 Z "i - 9 Z v 2 , vi - u 2 ), 

we have 

J s 

(P((fe-V)^)(r), <W r ,(r)) ro _ 1 , ro+1 dr + 2 / (P((v 2 • V)50y*)(r), 80 Y « (r)) m _i, m+1 dr 



+ IIMOILII^WII^dr + jf Ib2(r)||l||^(r)||4dr 

which together with the Gronwall-Bellman inequality, implies 



sup \\S9 Y ,( S )\\ 2 m + - \\Se z ,(r)\\ 2 m dr< C(u)R"e K 1 e|Nlc([r- E ,T];ff™)- 

Therefore, if we choose e to be small enough, the solution map 9 Y v is a contraction mapping on 

the complete metric space C/Jj and then through a bootstrap argument, we obtain a unique function u £ 
C loc ((T ,T];H™) nLl c (T ,T;H™+ 1 ) satisfying (6 Y u,0 z «) = («, Vu) on (T ,T] x M d with T depending 
on u,T, ||G||m and ||/||L 2 (o,T;ff m - 1 )- From Proposition 14.31 and the contraction mapping principle, we 
have 

(X,Y,Z,Y ) -.= (x a ,r Q ,z s ,? fi ) 

is the unique local iJ m -solution of FBSDS (jl.3p and all the other assertions hold as well. 

From Remarks 13.31 and 13 .41 we deduce that there exists some p £ Lf oc (To, T; H™) such that Yq = Vp. 
For each t £ (Tq,T] a.e. x £ R d , define the following equivalent probability Q l ' x : 

dQ^^expf--^ 9 Y {s,X s (t,x))dW s --^ £ \b( s ,X s (t,x))\ 2 ds^j dP. 

Then we have 

' dX s {t,x) = y/vdW' s , s€[t,T}; 
X t (t, x) = x; 

-dY s (t,x) = [(9 Y ■ V)9 Y +f + Vp] {s,X a (t,x))ds - ^V9 Y (s,X s (t,x))dW^ s £ [t,T]; 
Y T (t,x) = G(X T (t,x)), 
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where (W',Q t,x ) is a standard Brownian motion. For any £ £ C£°(R) (E>C^°(0), Ito's formula yields that 



C(s,X s (t,x))=((T,X T (t,x)) 
and thus, 

E q t,*[(9 Y , C)(t,x)]+vE 

— Eft*," 



(d r + -A)((r,X r (t,x))dr 



V((r,X r (t,x))dW^ 



(C(T,X T (t,a:)), G(z T (t,z))) 
+ J' ((-dsC-^AC, 8 Y >) + (C, (0 • V)0 + Vp + /»(«,*.(*, a:)) <fe 



Integrating both sides of the last equality with respect to x, we have 

i-T 



(C(«), e Y (t)) = (C(T), G) + J [- (d s {(s), 9 Y (s)) + (C(s), -A6 Y (s) + (9 Y ■ V)M*))o] da 

Hence, (8 Y ,p) is a strong solution to Navier-Stokes equation (|3.16p (see [101 EI])- Because of the re- 
versibility of the above procedure and the uniqueness of the if m -solution of FBSDS (11.3[) . we prove the 
uniqueness of the strong solution for Navier-Stokes equation (|3.16[) as well. 



Remark 4.5. Proposition 14.31 together with the contraction mapping principle serves to guarantee the 
existence and uniqueness of the local i? m -solution of FBSDS (|1.3[) . Furthermore, in a similar way to the 
above proof, we can prove that 8 Y of Proposition 14.31 is in fact the unique local strong solution of the 
following PDE 

d t u + | Au + ({b + au) ■ V)u + = 0, t < T; 
u(T) = i,, 

which is the well-known Burgers equation if a — 1 and 6 = 0. 



(4.26) 



5 Global results 



5.1 The case of small Reynolds numbers 

We work on the d-dimensional torus T d = M. d / (L x where L > is a fixed length scale. Denote by 
(H™' q (T d ; K rf ), || • || nir T<i) the M d -valued Sobolev space on T d , each element of which is divergence free. 
For g = 2, write (H^'(T d ; R d ), || • ||„. T<J ) for simplicity. Let G € i?™(T d ;R d ) of zero mean for m > d/2. 
Consider 



dX s (t, x 
X t (t,x 
-dY s (t,x 
Y T {t,x 

-dY s (t, x 



= Y s {t,x)ds + sfvdW s , se[t,T]; 
= x; 

= Y (s,X s (t,x))ds- V^Z s {t,x)dW s , se [t,T]; 
= G(X T (t,x)); 
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- dM. 



3 Y t ^(t,x + B s )[B^ 



B\ 



(Bi 



Bi ds 



s e (0, oo); 



Y oo (t,x) = 0, 



(5.1) 
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where the external forcing is not introduced for simplicity. In a similar way to Theorem 13. 31 FBSDS (15.11) 
admits a unique local ff m -solution on some interval (Tq,T]. Moreover, we have 

\\0Y(t)\\ 2 m . td + yJ \\9 z (s)\\ 2 m . Td ds 

= \\G\\Lt« + 2 / (P(6z(s)6 Y (s)), e Y ( s )) m _ hm+1 . Td ds 

< ll G llm;T<* + C J \\dY(s)\\ m ;r40z(8)\\ m -l.jd\\0Y(s)\\ m+l fx*ds 

< \\Gfm;t* +C [ ¥Y{ S )\\ m ^Pz{s)\\ 2 m . Jd ds, 

where we have used the fact that by Poincare inequality and the scaling properties, 

\\d Y (s,-)\L;T« <cl||vM8,-)IL ; t*. se[t,T], 

with 9 Y {s, •) being mean zero and the constant C being independent of L. Thus, 

IMtCjr- + /V - CL\\e Y { s )\\ m , v )\\e z { s )\\l VJd ds < \\G\\ 2 m . Jd . 

If we take the Reynolds number R := Lv~ 1 \\G\\ m .f d < C" 1 , then for this local solution (X, Y, Z, Yo) we 
always have 

l|0r(*)llm;T' i < l|G ! |lm;T< i ) * e (T(,,T]. 

Using bootstrap arguments, the local solution can be extended to be a global one. In summary, we have 

Theorem 5.1. Assume that G e H™(J d ;R d ) (m > d/2) is mean zero. FBSDS (|5.ip admits one and 
only one local H m -solution (X, Y,Z, Yq) on some time interval (Tq,T] with By € H™(T d ) being spacial 
mean zero. Moreover, there exists a positive constant Rq (= A as above) such that if the Reynolds number 
R < Rq, our local H m -solution can be extended to be a time global one and for this global H m -solution 
we have 

||M*)llro;T<* < \\G\\ m . Jd , for any t € [0,T]. 

5.2 The two-dimensional case 

Consider the two dimensional case. For simplicity, we assume / = and m > 3. Then under the 
assumptions of Theorem 13 . 31 let (X, Y,Z, Y ) be the local i? m -solution of FBSDS (|1.3I) on the time 
interval (Tq,T}. Define the vorticity field: 

# Y :=Curl 6 Y := d x J Y - d x . 2 9 Y , 

which is scalar-valued. Consider the following FBSDE: 

' dX s {t,x) = 9 Y (s,X s (t,x))ds + ^dW s ; 
< Xt(t,x)=x; _ 

dY s (t,x) = ^Z s dW s ; V ' 1 

, Y T (t, x) = (Curl G)(X T (t, x)), T < t < s < T, x G M 2 . 

By Proposition 14. 31 and Theorem 13.31 we have 

Y (t,x) = Y t (t,x) = £[(Curl G){X T {t,x))] . (5.3) 

In view of [301 page 117, Proposition 3.8], we have 

\\8z\\l-<F) ^ C 0- + ln+ HMa + ln+ lhMo)(l + ||tMU~(*»))> ( 5 - 4 ) 
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where ln + y :— V \ny. 

On the other hand, since V • Qy = 0, we have 



det (VX s {t,x)) = exp U (V • 6 Y )(r, X r (t, x)) dr J = 1, 

which implies 

\\#Y(t)\\ L « m < ||Curl G|| L9(R2) < C\\G\\ m , q e [2,oo]. 
Thus, in view of (|5.4|) . we have 

\\8z\\l~(k>) < C(\\G\\ m )(l + ln+ \\e Y \\ m ). (5-5) 
From the identity equation and the estimate (|2.4[) of Lemma 12.21 we have 

IIM*)llm + «' f \\dz{r)\\ 2 m dr = \\G\\ 2 m + 2 f {6 z 9 Y (r), Q Y {r)) m dr 

J S J S 

< \\G\\l l + C £ pzir^^^pyir^dr 

which by Gronwall inequality implies 

||M*)IU < C||G|| m exp (j* ||^(r)|| i0 c (R2) dr^j . 

In view of (|5.5p , we have 

ln+ ||^(s)|| m < C(\\G\\ m ,T) ( 1 + / ln+ ||0 r (r)|| dr) . 



(5-6) 



Gronwall inequality yields that 

sup \\0Y(s)\\ m <C(\\G\\ m ,T), Vte(T ,T]. (57) 

s e[t,T] y ' 

using a bootstrap argument, we can extend the local ff m -solution (X, Y,Z, Yq) of FBSDS (|1.3p into a 
global one. Therefore, we have 

Theorem 5.2. Let d = 2, m > 3 and G € iJ™. TTien our FBSDS (TQ|) wrf/i / = admits a unique 
H m -solution (X,Y,Z,Y ). 

6 Appendix 

6.1 Proof of Lemma 13.11 



It is sufficient for us to prove (|3.3|) with Z = 1, from which (|3.4j) follows by Fubini Theorem. 

First, taking a nonnegative function ip g C^ >0 (]R <i ;K), we consider the following trivial FBSDE: 

dX r (t,x) =b(r,X r (t,x))dr + ^dW r , T-s<t<r<s; 
X t (t,x) = x; 

r- (6-1) 

dY r (t,x) = y f vZ r (t,x)dW r , r g [t,s]; 
k K T (t,a;) = ¥'(-y.(t,a;)). 



2G 



In view of Lemma [4. 2 1 and the proof therein, FBSDE (16.11) is a particular case with <j> = therein, and 
moreover, the assertions of Lemma [4.21 still hold for (|6.ip . as Lemma [3.11 will never be involved in the 
proof of Lemma [4721 if = 0. Therefore, for almost all x € R d our FBSDE (|6.ip admits a unique solution 



(X.(t,x),Y.(t,x),Z.(t,x)) € S 2 (t,s;R d ) x S 2 (t,s;R d ) x L^(i,s;M d ), 
and for this solution (X, Y, Z), there hold 

By € C([f, *]; ff m ) n L 2 (i, s; £P l+1 ) 
^(r,X r (t,a:))= V j(X.(t ) a:))-v^ / 0z(r, JT T (t, a;)) dW T , a.s. 

J r 

and 

0z(t,a:) = V9y(t,x), Y r [t,x) = 9 Y (r,X r {t,x)), Z r (t, a:) = dz{r, X r (t, x)), a.s.,t < r < s. 
In an obvious way, we have 

Y r (t,x) = E[tp(X s {t,x))\& r ) > 0, a.s.re [Mi- 
Define the following equivalent probability measure 

dQ*'* =; exp f / 6(r, X r (i, x)) dW s - \v~ x / |6(r, JC r (t, x))| 2 ds J dP. 



In view of (|672|) . FBSDE ([671]) reads 

dX r (i, a;) = y/vdW^., * < r < s; 
X t (t, x) = x; 

-dY r (t,x) = Z r (t,x)b(r,X r (t,x))dr - \fuZ r {t, x) dW' r 

= (6- V)6» r (r,X r (t,a;))dr - ^Z r (t, x) dW' r , r G [t,s] 
r a (t,a:)=v»(X a (t,a;)), 

where (VF', Q 4 ' 31 ) is a standard Brownian motion. Therefore, 
9y(r, x) dx 

Eqt,* [9y(r, X r (t, x))] dx 

Eqt,* [(fi(X s (t,x))] dx + 



IR d Jr 

ip(x)dx+ j I (b ■ Y)9 Y (t,x) drd.i 

tp(x) dx - 



Eqt, x [(b-W)9 Y (T,X T (t,x))} drdx 



< / <p(x) dx 



(div b)9y (t, x) drdx 
||div6(r)||ioo / 9y(r,x) dxdr, 



or 



> / tp(x)dx— I ||div6(r)||ioo / 9y(T,x)dxdT. 

JR d Jr 

Using Gronwall inequality, we have 

I 6 Y {r,x)dx< k- 1 I ^x)dx, We M 



(6.2) 



(6.3) 
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with 

K := e _l|divf,|ll - 1 ( . T i i00 ). 

Taking r — t, we have (|3.3[) for K := k^ 1 . 

For the general function <p g C^°(R d ;R) without the nonnegative assumption, we choose a positive 
Schwartz function h and a nonnegative function <p £ (R d ;R) such that 

supp<^ C tp(x) = 1} . 

Set 

ip £ := \J (p 2 + ship, for e e (0, 1). 
Then in view of the above arguments, we have 

< k\\v> e \\li(r*) < / E[\ip £ (X s (t,x))\]dx < K\\ip £ \\ L i {Rd) . 

Letting e — > 0, we conclude from Lebesgue dominant convergence theorem that (|3.3[) holds for all ^ € 
C£°(R d ;R). 

Finally, for any tp £ L 1 (R), we choose a sequence {y™, n G Z + } C C£°(R ; R) such that lim n _ j . oc ||y> — 
¥' T *IU 1 (R) = 0. Then, by (|3.3[) . {<^ n (X s (t, a;))} is a Cauchy sequence in L 1 (fi x R d ;R). It remains to show 
that (p(X s (t, •)) is the limit. 

Through the above approximation, we can check that (|3.3p holds for any continuous function of a 
compact support. Therefore, if A C R d is a measurable, bounded subset of zero Lebesgue measure, then 
the dP x dec-measure of the set {(w,x) € il x R d : X s (t, x) S A} is zero. Thus, the almost everywhere 
convergence of <p n to ip in R d implies that of ip n (X s (t, •)) to ip(X s (t, •)). 

Hence, cp n (X s (t, x)) converges to tp(X s (t,x)) in L 1 (f2 x R d ;dP x dx). Since (|3.3p holds for each 93", 
passing to the limit, (|3.3I) holds for any y> € L 1 (R). We complete the proof. 

6.2 Proof of Lemma 14.11 

Step 1. For the trivial case of a = 0, the proof is referred to [36j [37] and it becomes a special one in 
[U EOl [29l |35] . It is sufficient to consider the case of a 7^ 0. 

Step 2. We prove the existence of the solution (X, Y, Z) which satisfies all the assertions of Lemma 

Choosing two positive real numbers e < T and M whose values are to be determined later, define 
V£:= {« G C([T-e,T];iOnL 2 (T-e,T;tf™ +1 ): 

IMIc([T-e,T];ff™) + 3 " Vu ^(T~e,T;H^ + 1 ) - Af2 }- 

For each £ € V^-, through a similar and simpler approximating method to be used in Proposition 14.31 we 
prove that the following FBSDE: 

dX s (t,x) = [b(s,X s (t,x))+a((s,X s (t,x))] ds + ^dW s , T - e <t < s <T 
X t (t,x) = x; 

\ r- ( 6 - 4 ) 

-dY s (t,x) = <t>(s,X s (t,x))ds-y/vZ s (t,x)dW s , s£ [t,T] 
k Y T (t,x)=iP(X T (t,x)) 
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admits a unique solution triple (X^,Y^, Z 1 ^) which satisfies 
WWL + v f \\0z,{r)\\ 2 m dr 

J S 

= Uf m + [ 2{({b + aQ-V)0 Yi {r),6 Yi {r)) m - 1 , m+1 dr+ I 2{</>{r), 9 YC (r)) m -i. m +i dr 



IMl7„ ! C\ j T \\b(r) + aC(r)\\l\\8 YC {r)\\ 2 m dr+ C \\4>(r)f m _ x dr 



(6.5) 



T {¥Yc{r)\\ 2 m + \\e zi {r)\\ 2 m )dr, 



where we have used Remark l2.ll By Gronwall inequality, we have 



2 



sup ||^ C ( S )|^ + -/ \\6 zi {r)\\ z m dr 

sG[T-e,T] z JT-s 



(6.6) 



Choosing M to be big enough and e to be small enough, we obtain 

r-T 



sup \\0 YC ( s )\\ 2 m + !i f \\e z( (r)\\ 2 m dr<M 2 . 



On the other hand, for any £i, C2 G Vj^, setting 

(S9 Y( , 59 z c , 5Q := (9 Y d - 9 Y c 2 , @ z <i - z c 2 , Ci - C2), 

we have 

\\69 Y< (s)\\ 2 m + v£\\6e z c(r)\\ 2 m d r 
= 2 / a((6( ■ V)6 Yil (r), tffly-c (r)) OT _i im +i dr + 2 / a((C 2 ■ V)<S0 y c (r), £0 y <(r)) m _i, m+1 dr 



, /' 



S J / i|l«»i«(r)li;, + ll«««(<-)E,)* 

2 11/3 /„M|2 , , / II/- /_M|2 Mr/) / \||2 



C / ||*C(r)OWr)|&dr+ / ||C 3 (r)||4,||^c(r)||^dr 



which together with the Gronwall-Bellman inequality, implies 



V 



sup \\66 Yi {s)r m + - \\58 z t(r)\\ 2 m dr< CM'e M 1 s\\8C,\\h {[ T-e,T];H™y 

sE[T~e,T] 2 J T _ e 

Therefore, if we choose e to be small enough, the solution map $ : C t— ► is a contraction mapping 
on the complete metric space Vj^ and then through bootstrap arguments, we obtain a unique function 
C G C 1oc ((t,T];H™) n tf oc ( T) T;H?+ 1 ) satisfying (0 yf ,0 z f) = (C.VC) on (r,T] x R d . where T - r 
continuously depends on ||0||i,2(o,T;fl" m - 1 )! ll^l|c([o,T];if">), HV'llmj v and a. Furthermore, By Lemma B~2l 
solving FBSDE (f6~?]) with C replaced by C, we get (X^, , Z?) which is a solution to our FBSDE (gTTJ 
and satisfies (J4T2]), (gTJJ and (|44|) . 

Step 3. In a similar way to Step 2 of the proof of Proposition 14.31 we prove the uniqueness. We 
complete the proof. 
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